Evidence for decay of spin-waves above the pseudogap in underdoped YBa 2 Cu 3 06.35 
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The magnetic spectrum at high-energies in heavily underdoped YBa2Cu30e.35 (T e =18 K) has 
been determined throughout the Brillouin zone. At low-energy the scattering forms a cone of spin 
excitations emanating from the antiferromagnetic (0.5, 0.5) wave vector with an acoustic velocity 
similar to that of insulating cuprates. At high energy transfers, below the maximum energy of 
270 meV at (0.5, 0), we observe zone boundary dispersion much larger and spectral weight loss 
more extensive than in insulating antiferromagnets. Moreover we report phenomena not found in 
insulators, an overall lowering of the zone-boundary energies and a large damping of ~ 100 meV 
of the spin excitations at high-energies. The energy above which the damping occurs coincides 
approximately with the gap determined from transport measurements. We propose that as the 
energy is raised the spin excitations encounter an extra channel of decay into particle-hole pairs of 
a continuum that we associate with the pseudogap. 

PACS numbers: 74.72.-h, 75.25,+z, 75.40.Gb 



There is a close relationship in high temperature su- 
perconductors between the carrier density, the magnetic 
ordering and excitations and the unusual ci-wave pair- 
ingiii^ii 3 -^ The superconducting phase first appears when 
the concentration of carriers is about 5.5 % and this small 
concentration of defects is also able to destroy the long 
range antiferromagnetic order. One of the most sur- 
prising results in the underdoped region of the phase 
diagram is the occurrence of a marked decrease in the 
charge scattering rate of the optical conductivity below 
the surmised pseudogap energy! 5 - The pseudogap energy 
decreases with increasing doping and its origin is contro- 
versial with many models and theories proposed^^^ 
Understanding underdoped high temperature supercon- 
ductors clearly requires an understanding of this pseu- 
dogap and more generally of how the antiferromagnetic 
Mott insulator becomes a superconductor for such small 
doping. 

YBa 2 Cu3 06.35 (YBCOe+z with x=0.35) is a supercon- 
ductor (T c =18 K) with a planar hole doping of about 6%. 
Low energy neutron scattering has shown that at this 
doping there is strong quasi-elastic scattering, showing 
the existence of commensurate antiferromagnetic short 
range order, and that there is no well-defined resonant 
mode* 1 - 1 - In contrast in YBCO6.5, and for similar hole 
doping, with T c = 59 K and a doping of 9 %, the com- 
mensurate quasi-elastic peak is replaced by incommen- 
surate excitations for energies less than a resonance at 
33 meV where the excitations become commensurate at 
the antiferromagnetic zone centre (0.5, 0.5) in units of 
27r/aj 12 ' 13 i 14 ' 15 i 16 : 17 Since both these materials are super- 
conducting it follows that none of these low energy prop- 



erties of the magnetism are essential for the occurrence of 
superconductivity. Our paper represents the first obser- 
vation by a direct spectroscopic probe of the destruction 
of spins by the pseudogap states, and has revealed an 
energy range where this onset occurs. Based on a com- 
parison to cuprates in the overdoped region of the phase 
diagram, we propose that this decay of spin excitations 
is a common feature in all cuprate high-temperature su- 
perconductors. 

To date there are very few measurements of the spin ex- 
citations at energies comparable with the zone boundary 
excitations at about 2J in the parent antiferromagnetic 
insulators and none exist for the heavily underdoped su- 
perconducting materialsj 18 ' 19 i 20 We report neutron scat- 
tering measurements in this region for YBC06.35- We 
compare the results with those obtained for insulating 
two-dimensional antiferromagnets and show that the low 
energy magnetic excitations are very similar. The high- 
energy excitations are very different and exhibit strong 
damping and reduced spectral weight near the magnetic 
zone boundary. These properties are characteristic of 
the dispersion curve entering a continuum of excitations 
as found in magnetic metals and low dimensional mate- 
rials. We suggest that the contiuum is associated with 
the pseudogap states and the large damping occurs when 
the energy of the magnetic excitations is larger than the 
pseudogap. 

The sample consisted of seven twinned crystals with a 
total volume of ~ 7 cm 3 aligned with a combined rock- 
ing curve width of ~ 2° and such that Bragg reflections 
(HHL) lay within the horizontal plane. The orthorhom- 
bic lattice parameters of the crystals were o= 3.843 A 
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FIG. 1: Constant energy slices are shown. The wavevector 
transfer is given by (H, K) = (£ + rj, £ — rf). 
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&=3.871 A and c=11.788 A. The measurements were per- 
formed at high energies using the MAPS time-of-fhght 
spectrometer at the ISIS facility. A Fermi chopper was 
spun at 400 Hz and monochromated the incident beam 
with energies of 150 meV, 310 meV, and 450 meV. The 
effective resolution was selected by combining pixels and 
time channels using the MSLICE programmed Typi- 
cally the energy integration was chosen to be 10-30 meV 
(full-width at half maximum) while the integration in 
wave- vector transfer was between 0.025 and 0.1 recipro- 
cal lattice units in the two-dimensional plane. 

In Fig. [1] we show constant energy slices through the 
intensity distribution of magnetic neutron scattering for 
incident energies of 310 meV and 450 meV. The coloured 
contour plots show how the scattering depends on wave- 
vector transfer in the two-dimensional CuC>2 plane. In 
Fig. [TJ the incident neutron beam was parallel to the 
c-axis and the results have been symmetrised to take ac- 
count of the intensity arising at equivalent wave-vector 
transfers. At low energies, the scattering is confined to a 
small circular disc centred on the antiferromagnetic Bril- 
louin zone centre, (0.5, 0.5). With increasing energy the 
size of the circle increases and the centre is depleted so 
as to form a ring whose diameter increases linearly with 
energy. Above ~ 200 meV the pattern becomes more 
complicated with the scattering intensity displaced to- 
wards the zone boundaries. In Fig. [2] we show a series of 
cuts through the data either to give constant energy scans 
or to give constant wave-vector scans through the data 
so as to obtain the dispersion from the antiferromagnetic 



FIG. 2: The left hand panels show cuts at constant energy 
transfer and with Ei=310 meV. Constant-Q cuts taken with 
Ei=450 meV are shown in the right hand panels for momen- 
tum positions around the zone boundary. The data was in- 
tegrated over (AH, AK) = (±0.05, ±0.05) and a background 
was taken at (0.5, 0.5). 

peak position at (0.5, 0.5) as shown in FigJU 

To extract the energy dependence of the acoustic and 
optic spectral weight, we varied the angle ip, defined as 
the angle between the incident beam and the [001] axis, 
with the rotation axis along [110]. With the incident en- 
ergy fixed at 150 meV and angles of ^=0, 20, and 30 the 
variation in the intensity with L enables the optic and 
acoustic components to be separated by Fourier analysis 
knowing the bilayer spacing <i/c=0.28. This is the same 
analysis and method applied previously to the YBCO6.5 
superconductors 16 As shown in Fig. [3]c) At low energies 
below 40 meV, all the scattering occurs in the acous- 
tic channel. At higher energies the optic mode intensity 
increases and at 75 meV it is larger than the acoustic 
mode intensity as shown by the reversal of the maxima 
and minima compared with the results at E=35 meV. At 
larger energies the dispersion energies and intensities of 
both optic and acoustic modes are within error the same 
and similar to the insulating antiferromagnet. 

For energy transfers below 200 meV, the momentum 
integral of the scattered intensity is shown in Fig. [3] a) 
as a function of the two-dimensional wave vector. Linear 
spin wave theory with a constant renormalization factor, 
Z x =0.51 (Ref. \2w . gives the dashed line. The intensity 
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FIG. 4: The momentum disp ersion is compared with the in- 
sulator YBCOs.is (Ref. GJ) and YBGOe.s (Ref. [H). C is 
measured from the zone centre (0.5, 0.5). 



FIG. 3: a) The momentum integrated intensity is plotted as a 
function of energy transfer, b) The energy integrated intensity 
as a function of wave vector relative to (0.5, 0.5) around the 
zone boundaries, c) The intensity as a function of L used to 
extract the acoustic and optic intensities. 



loss at high energies while subject to uncertainties in the 
absolute normalization and in the estimation of the back- 
ground is likely a forerunner of the large zone boundary 
loss discussed below. For energy transfers below ~ 200 
meV, the magnetic excitation energies and intensities re- 
semble the insulating antiferromagnet YBCO6.15 and the 
acoustic and optic intensities are similar to those of the 
insulator. In particular, the slope of the dispersion curve 
gives a spin- wave velocity of 592 ± 30 meV • A, which is in 
agreement with the velocity of ~ 600 meV • A measured 
for insulating YBCOg.is*^ There are, however, consider- 
able differences from spin wave behaviour at high-energy 
transfers above the range of hco ~ 150-200 meV. 

Firstly, the energy widths of the excitations, shown 
in the constant Q plots in Fig. [2] are about 100 meV 
which are much larger than the experimental resolution 
(horizontal bar) of about 10 meV. Therefore, the zone 
boundary excitations have acquired strong damping. 

Secondly, we observe an overall lowering of the zone 
boundary energies that is inconsistent with that of a 
simple antiferromagnet with a coupling between local 
spins. The dashed line in Fig. 2] shows that a near- 
est neighbour model J=130 meV can fit the velocity at 
low-energies but fails to describe the dispersion at high 
energies above ~ 150-200 meV. The dotted curve shows 
that the La2Cu04 nearest and next-nearest interactions 
of 112 meV and -11 meV respectively, also overestimates 
the frequencies of the zone boundary excitations. The 



solid line represents an attempt to reproduce the zone 
boundary energies (with nearest and next nearest neigh- 
bour interactions of 96 and -10 meV respectively) and 
it shows that the data cannot be described by a consis- 
tent description with only near neighbour interactions. 
All calculations include a rcnormalization factor calcu- 
lated to be Z c =1.18£2- We also observe anisotropic be- 
haviour at the zone boundary almost twice as large as 
that found in the insulator La2Cu04. The zone bound- 
ary dispersion in YBCO6.35 carries the same sign as that 
observed in La2Cu04 and opposite to that in the layered 
Sr 2 Cu 3 4 Cl 2 and CFTD^i^ 

A third difference between our results and linear spin 
wave theory is shown in Fig. [3] b) which plots the inte- 
gral over all energies taken from constant-Q cuts above 
150 meV. The zone boundary intensity is about a fac- 
tor of two less than predicted by linear spin wave theory 
with a renormalization factor (Z x =0.51) i 19 i 26 Figure [3] 
b) also shows that the difference decreases on approach- 
ing the magnetic zone centre. A decrease in the zone 
boundary magnetic neutron scattering intensity has also 
been suggested based on measurements in superconduct- 
ing Lai.86Sro.i4Cu04i^ We note that in the square lat- 
tice CFTD, the spectral weight loss occurs only near the 
point (0.5, 0) whereas we observe a loss of spectral weight 
across a broader range in momentum^ 

The large energy widths, reduced zone boundary en- 
ergy, and suppressed spectral weight suggest that the 
magnetic excitations interact with other high energy ex- 
citations that do not scatter neutrons strongly. The mag- 
netic excitations could decay into two spin-waves but we 
find that it is not possible to conserve wave vector and 
energy with the creation of pairs of uncoupled magnetic 
excitations having the dispersion of Fig. @J2£ We note, 
however, that if the single-magnon linewidth is increased, 
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then the kinematic constraints are relaxed and the decay 
into a pair of magnons may be possible. It is also possible 
that in a S=l/2 system, there are two spin-wave states 
that are coupled and that lie below the two spin wave 
continuum. The interaction between these states and 
the magnetic excitations might give rise to the damping. 

The fact that the unusual high-energy dispersion can- 
not be made consistent with the spin-wave models that 
describe the low-energy velocity and intensity (and also 
the entire spin-wave dispersion in the insulators) points 
to the need for an energy dependent decay process. It is 
likely that the damping arises from the decay into elec- 
tronic particlc-holc pairs since the most energetic spin 
excitations lie above estimates for the pseudogap. One 
suggestion from thermal conductivity measurements is 
that the quasi-particle gap for heavily underdoped YBCO 
is ~ 160 meV, which is comparable to the energy where 
the magnetic excitations show strong damping^ 7 - The en- 
ergy is also comparable to the expected pseudogap energy 
measured by NMR and optical absorption^ Hence, for 
energy transfers greater than the pseudogap, the decay 
of a magnon into a particle-hole pair is allowed. 

In ferromagnetic metallic systems characterized by a 
continuum of particle- hole excitations (for example in Fe, 
Ni, or MnSi), spin- waves become strongly damped when 
both q and u> match those of the particle hole pairs i 31 i 32 
This has been suggested to account for the anomalous de- 
crease in intensity in YBCOg.85 for energies ~ 55 meV^ 
A similar decay process may account for strong damping 
of the high-energy excitations in YBCO6.35. It therefore 
appears that the large damping of the spin-excitations at 
high-energies maybe a common trait amongst all high- 
temperature superconductors across the entire phase di- 
agram unlike the resonance and incommensurate scatter- 



ing which have not been observed in all superconductors. 
It would be interesting to model this effect however, the 
t-J model has been applied to heavily doped supercon- 
ductors, but to our knowledge predictions for the lightly 
doped superconductor have not been made, possibly be- 
cause of uncertainty in the Fermi surface . 33 i 34 i 35 i 36 i 37 i 38 

Near the onset of superconductivity, the low energy 
spin response is strong and gives the appearance of spin- 
waves while at high energy we have shown that the spin 
excitations enter a region where the decay, decrease in 
energy, and intensity do not fit into a spin-wave descrip- 
tion. The evidence points to the presence of a particle- 
hole continuum as a possible explanation of the results. 
On this basis, we consider that the pseudogap is identi- 
fied to be in the range 150-200 meV for £=0.35, T c =18 
K. The damping observed for x=0.85, T c =85 K also sug- 
gests that the spin excitations decay near 55 meV and 
if this has the same origin it implies that the pseudogap 
has decreased with increasing doping. In contrast, the 
low energy results between £=0.35 and £=0.5 (Fig. |4|) 
show the presence or absence of quasi-elastic scattering, 
commensurate or incommensurate excitations, and the 
presence or absence of a resonance suggesting that these 
features are not present in all superconductors. We there- 
fore suggest that the existence of the high energy con- 
tinuum seems to be the feature which is common to all 
superconductors at least up to optimal doping, whereas 
the low energy features seem less universal. 
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